Ab initio study of semiconducting carbon nanotubes adsorbed on the Si(lOO) surface: 
diameter- and registration-dependent atomic configurations and electronic properties 
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We present the first ab initio study of semiconducting carbon nanotubes adsorbed on the un- 
passivated Si(lOO) surface. Despite the usual gap underestimation in density functional theory, a 
dramatic reduction of the semiconducting gap for these hybrid systems as compared with the elec- 
tronic gaps of both their isolated constitutive components has been found. This is caused by the 
changes in the electronic structure as the surface reconstructs due to the tube's proximity, the con- 
comitant electronic charge transfer from the nanotubes, and the band hybridization with silicon and 
carbon states resulting in the appearance of states within the energy gap of the formerly isolated 
nanotube. Furthermore, it is determined that semiconducting nanotubes exhibit weaker adsorption 
energies and remain at a greater distance from the Si(lOO) surface as compared to metallic nan- 
otubes of similar diameter. This effect may be useful for the solid-state separation of metallic and 
semiconducting nanotubes. 

PACS numbers: 68.35.-p,68.43.-h,68.43.Bc,73.22.-f 



I. INTRODUCTION 



First-principles studies illustrating the effect of techno- 
logically relevant semiconductor surfaces such as InAs,(l|] 
GaAs|2| and Si(lOO) on the electronic properties of single- 
wall carbon nanotubes (SWNTs) have been published 
recently. In the case of Si(lOO), the study was focused 
on determining the lowest-energy structural configura- 
tion and modifications to the electronic structure of a 
metallic (6,6) SWNT as a result of its interaction with 
this surface. This calculation was performed for a nan- 
otube in the proximity to an either clean|3, 3 or selec- 
tively hydrogen-passivated[3 Si(lOO) surface where the 
nanotube axis was parallel to the trench between adja- 
cent Si dimer rows. Remarkably, there have been no 
studies of this system for semiconducting tubes, in dif- 
ferent geometrical configurations, nor on the dependence 
of the properties of the hybrid system against nanotube 
diameter. These are relevant issues to be addressed, as 
experimental techniques permitting the ultra-clean de- 
position of SWNTs onto doped Si(lOO) and other semi- 
conductor surfaces in ultra-high vacuum (UHV) at room 
temperature have been reported. [^0,13 As contaminant- 
free atomistic manipulation becomes more feasible, the 
promise of molecular systems with tunable electronic and 
mechanical properties becomes a reality. 

It is clear that near-term applications for carbon nan- 
otubes in electronic and opto-electronic devices would 
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benefit from their integration with conventional semi- 
conductor platforms such as Si or GaAs. In this direc- 
tion, rectifying carbon nanotube-silicon heterojunction 
arrays have been demonstrated, and semiconductor 
heterostructures of GaAs/AlAs and GaAs/MnAs have 
been used as electrical contacts to individual SWNTs. |0| 
Recent experiments!ll| showing the existence of two pref- 
erential directions for the growth of SWNTs on silicon 
surfaces indicate that there is a non-negligible interaction 
between the SWNT and its substrate. One relevant issue 
which remains to be addressed is the nature of the me- 
chanical and electronic properties of the semiconductor- 
nanotube interfaces. 

The subject of this work is the interaction between the 
unpassivated Si(lOO) surface and semiconducting carbon 
nanotubes using density functional theory (DFT). l3| 
Those studies are complemented with additional calcula- 
tions on metallic nanotubes of comparable diameter. We 
study nanotubes that are in parallel or perpendicular ori- 
entations relative to the Si dimer row direction. We have 
found a striking and somewhat counter-intuitive reduc- 
tion of the semiconducting gaps for these hybrid systems 
which arc composed of two semiconductors. The nan- 
otubes presented in this study have diameters between 5 
and 12 A. This work was comprehensive in order to find 
robust properties on these hybrid systems; in particular 
those independent of diameter and chirality as experi- 
ments do not have fine control of these variables yet. 

Section|n]discusses the theoretical approximations em- 
ployed and the structural configurations considered in 
this study. In Section ITTll we show the resulting atomic 
configurations of the combined Si(100)-SWNT system, 
which have marked trends depending on the original 
electronic character of the SWNTs involved (metallic or 
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TABLE I: Lattice constants obtained with our bases, to ex- 
emplify the flexibihty of our basis set. 



Element 


Theoretical estimate (A) 


This work (A) 


C (Diamond) 


3.570T 


3.486 


C (Graphene) 


a 2.450* 


a 2.448 




c 6.500* 


c 6.516 


Si (Diamond) 


5.430* 


5.401 



fReference |2U 
^Reference l2a 



TABLE II: Carbon nanotubes with length in their unit cell 
commensurate within 2% to the Si(lOO) surface's supercell. 
In bold, the respective values for the (6,6) SWNT previously 
studied by Orellana et al. 3, 4, 5] 

Nanotube Diameter (A) Length (A) DFT Gap (eV) 



(6,2) 


5.657 


15.391 


0.860 


(8,4) 


8.302 


11.294 


0.945 


(12,4) 


11.314 


15.391 


0.699 


(9,3) 


8.486 


15.391 


0.000 


(6,6) 


8.153 


2.465 


0.000 



semiconducting). Results from Voronoi and Hirshfield 
population analysis _^13j are also provided and show elec- 
tronic charge being transferred from the nanotubes to 
the silicon slab. The resulting band structures at equi- 
librium, as well as the projected density of states (PDOS) 
over the atomic species involved are reported also. The 
PDOS indicates a clear hybridization of the bands be- 
tween surface and nanotube states. This hybridization 
is further confirmed by depicting wavefunctions near the 
electronic gaps. Finally, we present the change in energy 
prior to the force minimization procedure, as the nan- 
otubes are rotated about their axes in close proximity to 
the Si(lOO) surface. Conclusions are presented in Section 



II. METHODS 

Our calculations in the local density approximation^^ 
(LDA) were performed with the SIESTA code." 
The exchange-correlation potential employed is the 
one parameterized by Perdew-Zunger,|16j| based on the 
Ceperley- Alder data.0| Core electrons are replaced by 
norm-conserving TrouUier-Martins pseudoDotentials.|l8j 
For greater variational freedom, a double-^ basis set for 
s and p orbitals, and a single-C basis set for d orbitals 
was constructed using the prescription of Junquera et 
al. 19] To ensure the flexibility in our basis sets, their 
parameters were optimized by means of the simplex al- 
gorithm (20| on graphite and diamond for the carbon 
basis, while the silicon basis was optimized only in the 
diamond structure. As the hydrogen atoms only served 
to passivate dangling bonds, the hydrogen basis was not 
optimized. As shown in Tabled the lattice constants ob- 
tained from the optimized bases compare well with pre- 
vious theoretical estimates. We use the p(2x2) recon- 
struction of the Si(lOO) surface for calculations involving 
semiconducting SWNTs. The gap obtained is of 0.235 
eV. In all cases a slab with six silicon monolayers with a 
height of 7.81 A was employed. The bottommost layer 
is hydrogen(II)-passivated forming a dihydride arrange- 
ment. The SWNTs will be placed in proximity of the 
uppermost unpassivated layer. The area of the unit cell 
is L X L with L = 7.639 A. The vacuum region in the 



vertical direction is fixed in order to provide at least 10 
A separation between periodic images when the SWNTs 
are in place. 

Since we are interested in trends that could comple- 
ment both experimental work in this area and published 
theoretical results on metallic SWNTs on this surface the 
SWNTs chosen satisfy the following criteria: 

1. They are semiconducting. 

2. Their diameter is of the order of 10 A. 

3. They are commensurate with the underlying sur- 
face. 

While (n,n) (metallic) SWNTs happen to be commensu- 
rate to within 3% with the underlying Si(lOO) surface, Q 
the shortest commensurate semiconducting SWNT, of in- 
dexes (3n,n) -and 3n — n ^ 3q with q integer- would 
be 15.391 A long. The large unit cell sizes involved 
in calculations have deterred researchers from perform- 
ing calculations in this system with semiconducting nan- 
otubes. Table |nl lists the nanotubes studied, their di- 
ameters, lengths as well as their DFT gaps. Tubes with 
chiral indexes (2n,n) are 11.294 A long, so a supercell 
constructed out of two nanotube unit cells is required in 
this case to meet the surface's supercell length. Thus the 
slab supercells required for calculations involving both 
the surface and semiconducting nanotubes have surface 
areas equal to 3L x 2L or 3L x 3L, depending on the 
length of the SWNT, its diameter and the surface re- 
construction employed. This choice leaves more than 10 
A between nanotube images. In order to test the accu- 
racy of our calculations against published results, 0, Q 
benchmark calculations on a chiral metallic SWNT in- 
teracting with c(4x2) reconstructed Si(lOO) surface were 
performed. The surface gap in this case is equal to 0.223 
cV. The c(4x2) reconstruction minimizes the total en- 
ergy the most and therefore is the best candidate for the 
surface reconstruction. |2^ The metallic (9,3) SWNT is 
comparable in diameter to the (6,6) SWNT studied ex- 
tensively by Orellana et al. |^, ^ Due to the surface re- 
construction used and the requirement for the nanotube 
images to be separated by at least 10 A this required an 
even larger, AL x 2L area for the underlying surface. This 
benchmark calculation is also used to assess the effect 
of chirality and the exchange-correlation approximation; 
LDA in this work and GGA in Refs. H H |H 
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TABLE III; Averaged displacement (over the unit cell) in the 
z-direction (A) of silicon atoms closest to the nanotube; c.f. 
Fig-0 (^i) ~ i^i) refer to the parallel configuration and (/15) 
to the tubes in the perpendicular configuration. 



Nanotube: 


(6,2) 


(8,4) 


(12,4) 


(9,3) 


(hi) 


0.15 


0.11 


0.03 


0.14 


(h2) 


-0.18 


-0.29 


-0.41 


-0.44 


(h3) 


0.20 


0.15 


0.09 


0.43 


(h4) 


-0.05 


-0.06 


-0.04 


0.00 


(hs) 


-0.13 


-0.26 


-0.16 


-0.36 



III. RESULTS 

Figurenshows hybrid structures with maximum resid- 
ual forces less than 0.02 eV/A. The relevant geometric 
parameters are the minimum distance D between the 
SWNT and the dimer trench, the vertical atomic dis- 
placements (hi) resulting from the reconfiguration of the 
silicon atoms closest to the SWNT, and the azimuthal ro- 
tation angle (p. We want to emphasize the shorter bond 
lengths between the carbon and silicon atoms when the 
metallic tubes are involved, as well as a marked distor- 
tion of the nanotube, emphasized by the inner (red) cir- 
cle. In comparison, semiconducting tubes are in equilib- 
rium farther away from the surface and do not change 
their circular profile, as evidenced by the (green) inner 
circle perfectly matching the nanotubes' perimeters. To 
contrast with published work,[30| due to chirality the 
metallic nanotube in either configuration studied showed 
only seven bonds of varying lengths (2.03-2.27 A in the 
parallel configuration and 2.02-2.15 A in the perpendic- 
ular configuration). For the semiconducting nanotubes, 
the shortest distance between carbon and silicon atoms 
turned out to be 2.62 A and 2.87 A for the (6,2) tube 
in the parallel and perpendicular configuration, respec- 
tively. All the semiconducting tubes with larger diame- 
ters had their carbon atoms more than 2.84 A away from 
the closest surface silicon atom and in most cases more 
than 1 A above the distance in the C-Si bond of silicon 
carbide. This is indicative of a weaker bonding for semi- 
conducting nanotubes adsorbed on the Si(lOO) surface. 
For chiral nanotubes, the optimal geometrical configura- 
tion can not be known a priori, but it has to come out 
after a sweeping through all angles 0; this angular depen- 
dence is not as pronounced for semiconducting SWNTs. 
A detailed discussion follows. 



A. Optimized atomic configurations 

An equivalent plane-wave cutoff of 200 Ry was used 
to calculate the charge density on the real-space grid. 
To obtain the optimal configuration of these systems, we 
used a grid with a single fc-point (the F-point), except for 
the (6,2) nanotube, where a 2x2x1 Monkhorst-Packp^ 
(MP) grid was already employed for obtaining the relaxed 



structures. The adsorption energies were then calculated 
with a 2x2x1 MP grid. 



1. Distance from the SWNTs to the silicon surface and 
surface reconfiguration 

Total energies for the hybrid SWNT-Si(lOO) system 
with the SWNTs listed in Table |n] were calculated as 
a function of the distance from the nanotubes to the 
bottom of the the dimer trench. The optimal distance 
D is reached when the total energy is at the global 
minimum. While holding the bottom hydrogen layer 
fixed, the hybrid structures were relaxed by a conjugate- 
gradient method to minimize the residual forces until 
they no longer exceeded 0.02 eV/A. We also calculated 
the optimal distance D for (8,8) and (10,10) SWNTs for a 
fixed angle (j) in order to visualize and distinguish trends 
between semiconducting and metallic SWNTs. No re- 
laxation procedure was performed on the latter tubes. 
Under relaxation those distances would decrease slightly, 
but the trend indicated here will stand. It should also be 
mentioned that the trend obtained in Fig. [21 for metallic 
tubes remained the same when the c(4x2) reconstructed 
surface was replaced by a p(2x2) reconstructed surface. 
This means that the distinct trends observed in Fig.Qlare 
due to the different electronic character of the nanotubes 
involved, and not to the surface reconstruction employed. 
This is stressed in Fig. 13 where the structural configura- 
tion of the (9,3) tube and the p(2x2) reconstructed sur- 
face with maximum residual forces down to 0.02 eV/A, 
appears extremely similar to the one found for this nan- 
otube on the c(4x2) reconstructed surface. The diameter 
of the nanotube as well as its electronic character (metal- 
lic, semiconducting) are the determining factors on the 
relative distance between the nanotube and the surface. 
As the nanotube diameter increases, the distance D from 
the SWNTs to the bottom of the dimer trench also in- 
creases asymptotically towards 5.52 A, which is the lim- 
iting distance from the bottom of the Si dimer trench 
to a graphene plane (nanotube diameter — > 00), indi- 
cated by the gray line in Fig. |21 30] The results in Fig. [21 
show that metallic nanotubes exhibit a stronger diame- 
ter dependence on their optimal distance to the surface's 
slab, ml In previous work, the same distance to the sur- 
face was used for the (6,6) SWNT and a tube of 100 
A in diameter, in an attempt to explain experimental re- 
sults of multiwalled SWNTs on H-passivated surfaces, psj 
We speculate that the adsorption energy of nanotubes 
on an H-passivated surface should be lower than the one 
corresponding to an unpassivated surface, because the 
dangling bonds of the silicon surface are hydrogen passi- 
vatcd. This effect is worthy of further investigation, but 
is not included in this present study. For the semicon- 
ducting nanotubes studied the diameter dependence is 
not as marked. Yet, Fig. [21 shows that a semiconducting 
SWNT in the perpendicular configuration will be about 
0.5 A farther away from the silicon surface in comparison 
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FIG. 1: (Color online) SWNTs on the Si(lOO) surface in two configurations: parallel over the dimer trench, (a)-(b) and 
perpendicular to the Si(lOO) surface dimer trench, (c)-(d). D is the distance from the bottom of the SWNT to the bottom of 
the Si(lOO) surface dimer trench obtained in the total energy minimization. The vertical displacements (hi) involved in the 
surface reconfiguration and 0, the angle of rotation of the SWNT about its axis, are also indicated. An evident 'squashing' 
effect -an elongation along the z-direction- can be observed for the (9,3) nanotube but it was not present in any semiconducting 
nanotube. Note that for chiral tubes, a rotation about their axis and a displacement Al (not shown) along the nanotube's axis, 
are linearly dependent. To minimize the cell size, a p(2x2) surface reconstruction was employed when placing semiconducting 
nanotubes, and the c(4x2) surface reconstruction for the (9,3) nanotube. 



with its equilibrium distance in the parallel configuration, 
in accordance with experimental observation. |26j 

The electronic properties of the Si(lOO) surface are 
closely related to its atomic configuration. Table lTTll gives 
the average vertical displacements (hi) that the silicon 
atoms closer to the SWNT undergo due to the proxim- 
ity of the SWNT. The Si(lOO) surface in the proximity 
of the SWNT is subject to a striking atomic reconfigu- 
ration. For a given SWNT the amount of surface recon- 
figuration is always greater in the parallel case. This is 
consistent with the fact that tubes in the parallel config- 
uration stay closer to the surface (Fig. ^ and therefore 
interact more strongly than in the perpendicular case. 
We have to stress here that after atomic relaxation the 
semiconducting nanotubes undergo almost no structural 
change but remain very stiff. This is not the case for the 
metallic nanotube, as previously indicated. 



2. Relative strength of the C-Si bonds 

With the goal of visually expressing the varying 
strengths of bonds formed between carbon and silicon 
atoms for the metallic and semiconducting SWNTs, we 
show in Fig. 0| charge density isosurfaces for the (9,3) 
and (8,4) SWNTs. Results for the remaining nanotubes 
can be found in Table IIVI Very dense bonds can be 
seen for the (9,3) SWNT in either configuration with the 
0.37 e/A^ isosurface. The charge density isosurfaces show 
that the C-Si bonds are as strong as the Si-Si and C- 
C bonds in the surface and the nanotube, respectively, 
as can be iirferred from the similar width. In order to 
understand the varying strengths of the bonds and the 
role of chirality on bond formation, refer to Fig. de- 
picting the (9,3) SWNT in the parallel configuration. In 
this figure, surface dimers are located in the (x,z) plane. 
For the shortest bond (2.03 A) the relative coordinates 
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FIG. 2: (Color online) Distance D from the bottom of the 
dimer trench to the lowest carbon atom as a function of nan- 
otube diameter for fixed (f). Semiconducting SWNTs parallel 
to the dimer trench will be closer to the underlying surface by 
about 0.5 A as compared to semiconducting SWNTs in the 
perpendicular configuration for this diameter range. Notice a 
more pronounced diameter dependence and an overall smaller 
distance to the surface for metallic SWNTs. Dashed lines are 
drawn as a guide to the eye to facilitate the visualization of 
trends. 
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FIG. 4: (Color online) Charge density isosurface plots to visu- 
alize the relative strengths of C-Si bonds. Values of the densi- 
ties associated with each isosurface are shown. The strongest 
bond occurs for the metallic SWNT. 




FIG. 3: (Color online) The properties of the hybrid system 
do not depend on the underlying surface reconstruction as 
evidenced by the (9,3) nanotube on the p(2x2) reconstructed 
surface. The distance between the nanotube and surface is 
very close to that found for this nanotube on the c(4x2) re- 
constructed surface, see Fig. |5| Notice also the vertical elon- 
gation appearing again and the overall structural similarity 
to the structure presented in Fig. 



TABLE IV: Charge density (e/A^) at the bonds between nan- 
otube carbon and surface silicon atoms. Notice the smaller 
density for all studied semiconducting SWNTs in comparison 
with the bond densities for the metallic nanotube. 



Nanotube: 


(6,2) 


(8,4) 


(12,4) 


(9,3) 


Parallel 


0.235 


0.175 


0.120 


>0.370 


Perpendicular 


0.155 


0.130 


0.115 


>0.370 



of the involved carbon and silicon atoms are (Ax, AZ, 
Az)=(-1.00, -0.23, 1.75) A. The second closest carbon 
atom to a surface silicon atom has relative coordinates 
equal to (0.99, 0.35, 1.76) A, giving a bond length of 
2.05 A. Although the cartesian projections |Aa;| and Az 
are very similar in both cases, chirality breaks the mirror 
symmetry with respect to a plane formed by the SWNT 
axis and the bottom of the dimer trench once all bonds 
are considered (recall that only seven -instead of eight- 
bonds were formed), but it also shifts the carbon atoms 
by varying distances A/, according to their positions in 
the SWNT unit cell. Contrast this to the more symmet- 
ric configurations previously studied for the (6,6) SWNT 
in which pairs of bonding atoms have the same AZ and 
a similar value for |Aa;| for all bonds due to mirror sym- 
metry. 

For the (8,4) SWNT we had to lower the density to 
0.175 e/A"^ to identify bond formation in the parallel 
configuration. Lower density through bonds will reflect 
the weaker adsorption energies found for semiconducting 
SWNTs. Notice that this isosurface does not form bonds 
in the perpendicular configuration. It is at the lower 
0.130 e/A^ density when a single isosurface connecting 
the nanotube and slab can be seen. This is about a 
third of the density present in the strongest C-Si bond for 
the metallic SWNT. In conclusion, bond formation will 
be weaker for semiconducting SWNTs, and the weakest 
bond occurs for semiconducting SWNTs in the perpen- 
dicular configuration. These results are consistent with 
all findings of previous subsections. 

3. Amount of charge transferred to the silicon slab 

The results from Voronoi and Hirshfield deformation 
density charge analysis [l^] are summarized in Table Ivl 
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TABLE V: Amount of electronic charge per unit length (e/A) 
transferred from the SWNT to the slab from the Voronoi (Hir- 
shfield) deformation density analysis. 



Nanotube 


Parallel 


Perpendicular 


(6,2) 


+0.052 (+0.046) 


+0.033 (+0.029) 


(8,4) 


+0.055 (+0.049) 


+0.035 (+0.033) 


(12,4) 


+0.072 (+0.065) 


+0.035 (+0.031) 


(9,3) 


+0.064 (+0.083) 


+0.058 (+0.043) 



Because of the way they are defined, the results from 
those methods are independent of the numerical orbital 
basis, and they are also more meaningful physically: The 
Voronoi deformation density charge for a given atomic 
nucleus, for instance, is the absolute charge defined in 
the volume defined by all the points closer to that atomic 
nucleus than to any other nuclei. In all instances, the 
charge is transferred from the nanotube to the slab, and 
the magnitude of the charge transferred is larger when 
the tubes are parallel and on top of the dimer trench. 



B. Ab initio electronic properties of the hybrid 
Si(100)-SWNT system 

1. Band structures and densities of states 

In order to have reliable band structures a stringent 
force relaxation was performed. As mentioned in Section 
nil Al the atomic positions were subject to a force mini- 
mization procedure with a single fc-point (the T point), 
with the exception of the (6,2) SWNT, where a 2x2x1 
MP grid was employed. The relaxed atomic positions 
from this calculation were then used in computing the 
band structures and PDOS. The band structures were 
obtained with a 4x4x1 MP grid, while the PDOS was 
computed from at least a 12x12x1 MP /c-point grid (for 
example, the results for the (12,4) nanotube in the par- 
allel configuration were obtained with a 24x12x1 MP 
grid), using the converged self-consistent charge density 
obtained from the calculation with a 4x4x1 MP grid. 
The band structures and PDOS are depicted in Fig. 
In all the plots shown in this figure the SWNT was placed 
parallel to the F — J' direction. For the PDOS a Gaussian 
smearing of 10 meV half-height width and 300 sampling 
points in the (—1, 1) eV energy interval were employed. 
The Gaussian smearing is in most cases smaller than the 
electronic gaps obtained from the band structures. |32| 
The PDOS is computed with a finer mesh, independently 
from the band structure calculation and indicates the 
contribution of each atomic species to the band struc- 
ture. The two columns in Fig. [51 correspond to the re- 
laxed hybrid system in either the parallel or perpendicu- 
lar configurations. The most salient feature from Fig. [SJ 
as the band structures indicate, is the dramatic reduction 
of the gap for the hybrid systems composed of semicon- 
ducting nanotubes (with gaps of at least 699 meV when 



TABLE VI: Electronic gaps (in eV) for the hybrid systems 
composed of semiconducting nanotubes on Si(lOO) surface. 
We also include the (9,3) tube for completeness. 



Nanotube: 


(6,2) 


(8,4) 


(12,4) 


(9,3) 


Parallel 


0.099 


0.008 


0.033 


0.000 


Perpendicular 


0.083 


0.013 


0.040 


0.000 



isolated) and the Si(lOO) surface (with a gap bigger than 
200 meV) . The biggest gap found occurred for the system 
involving the (6,2) nanotube, and is equal to 83 meV. Re- 
fer to Table IVll for a comprehensive list of the electronic 
gaps. The projected densities of states in Fig.Oindicate a 
contribution from the carbon orbitals to the band struc- 
ture at energy values where the gap of isolated tubes is 
expected, highlighted by the arrows on the carbon PDOS 
(note the different scales for each species in the PDOS). 
We believe this reduction of the gap for the hybrid system 
with respect to the gaps of its constituent subsystems will 
occur even when more accurate quasiparticle (e.g. GW) 
calculations and experimental data are in place. Notice 
that, despite the strong structural modification observed 
for the metallic tube (Figs, and^), this is insuffi- 
cient in this case to open a gap in either configuration; 
refer particularly to Fig. where we zoomed the band 
structure about the Fermi level in the parallel case, to 
better appreciate this fact. The reader might appreciate 
the existence of a single band above the Fermi energy for 
the system involving the (8,4) SWNT in the perpendic- 
ular configuration. The PDOS indicates a couple of dips 
and the density at those dips is not equal to zero. The 
discrepancy here only comes due to the insufficient num- 
ber of fc-points employed. This system has the largest 
number of atoms, and we are certain a calculation with 
a finer mesh will result in a flatter PDOS in the region 
of interest, but we did not pursue further calculations 
due to their expensive nature. In the same note, the 
computational cost involved in computing a PDOS for 
the zoomed bands, which requires a finer smearing and a 
correspondingly increased number of fc-points, prevents 
us from performing such calculation. Nevertheless we are 
certain of the accuracy of the band structure, and there- 
fore of the accuracy of the electronic gaps for the hybrid 
system with a (8,4) tube adsorbed. Experimentally, the 
gaps for the nanotubes could be obtained accurately on 
the inert H-passivated substrate, prior to hydrogen de- 
passivation; the surface gap can be measured by a host of 
experimental techniques. The dramatic reduction of the 
gap in this hybrid system as well as the distinct struc- 
tural trends for semiconducting and metallic nanotubes 
are the main findings of our investigations. 

In Ref. I23 a gap opening for the small-diameter and 
achiral (5,5) nanotube was found in the perpendicular 
configuration. Due to a lesser symmetry in chiral tubes, 
which translates in longer bond lengths and a more com- 
plex atomic reconstruction and band hybridization, this 
effect is not seen in the (9,3) nanotube reported on this 
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FIG. 5: (Color online) Ab initio band structures and projected densities of states for SWNTs in different alignments with 
respect to the Si(fOO) surface. The first column shows the band structures and PDOS when SWNTs are aligned on top and 
parallel to the dimer trench, while the second column depicts band structures when the nanotubes are perpendicular to the 
dimer trench. Those results are obtained after atomic relaxation was performed. The ab initio results show a drastic reduction 
of the semiconducting gap for the hybrid system composed of semiconducting tubes and the Si(lOO) surface and a high degree 
of band hybridization. The different positions for the J, J' points reflects the difference in size of the unit cells considered. 
The band structures involving the metallic nanotube (in blue) show no gap opening at the Fermi level. The band structures 
are calculated with a 4x4x1 MP grid, while the PDOS was obtained with at least a 12x12x1 MP grid. Arrows indicate the 
contribution to the band structure from carbon atoms in the vicinity of the Fermi level. 



paper, and certainly more systematic studies would be 
required to address such a rather interesting effect; par- 
ticularly with regards to knowing if it is robust enough 
to be observed in other nanotubes. Our work indicates 



that any semiconducting nanotube on this semiconduc- 
tor surface results in a hybrid system with a gap smaller 
from that of the isolated constituent subsystems. 
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FIG. 6: (Color online) Band structure for the (8,4) and (9,3) nanotubes around the Fermi energy to visualize gap openings or 
lack thereof. We omit zooming into band structures for which the gaps can be directly visualized from Fig. |^ 



2. Highest occupied, and lowest unoccupied electronic states 
for the Si(100)-SWNT system 

Wavefunctions for electronic states at the r-point that 
are just below and above the Fermi energy -the (F- 
point) Highest Occupied and Lowest Unoccupied Elec- 
tronic States (HOES/LUES)- are depicted in Fig. □ 
From those plots we observe that the electrons respon- 
sible for metallicity are either (a) extended over both 
the SWNT and the silicon slab; (b) predominantly over 
the silicon slab or (c) confined to the SWNT. We also 
found an instance -the (12,4) SWNT in the parallel 
configuration- in which the wavefunction is localized at 
the interface between the nanotube and the surface. In 
most cases involving the semiconducting nanotubes a 
more pronounced contribution from the nanotube to the 
HOES is seen in contrast with a more surface- like LUES. 
Those findings can be understood in terms of the relative 
contributions to the band structure from the carbon and 
silicon orbitals as seen from the PDOS in Fig. [5| 



C. Adsorption energies 

Adsorption energies on the Si(lOO) surface for the car- 
bon nanotubes in this study are presented in Table IVIll 
The adsorption energies are obtained as the difference be- 
tween the total energies in the relaxed hybrid structures 
and those of the fully relaxed tubes and surface in the 
same supercell separated by 10 A. Notice the constant 
adsorption energy for semiconducting nanotubes in the 
perpendicular configuration, and the decreasing adsorp- 
tion energy with increasing diameter for semiconducting 
tubes in the parallel configuration. These two trends 
meet for a tube with diameter in between the (8,4) and 
(12,4) nanotubes, after which a semiconducting tube in 
the perpendicular configuration has a stronger adsorption 
energy. This is stressed by the bold font in Table IVIll In 
order to determine the best angular 'starting' structure 



TABLE VII; Adsorption energies per unit length (eV/A) af- 
ter full relaxation has been achieved. Notice the lower adsorp- 
tion energies obtained for semiconducting nanotubes. For the 
(12,4) tube (in bold) the perpendicular configuration turns 
out to be the most favorable. 



Nanotube: 


(6,2) 


(8,4) 


(12,4) 


(9,3) 


Parallel 


0.20 


0.18 


0.11 


0.63 


Perpendicular 


0.15 


0.15 


0.15 


0.60 



a chiral nanotube has to be rotated about its axis. Chi- 
rality precludes a highly symmetric configuration, poten- 
tially with the lowest energy, from being easily visualized. 
The energies as a function of 4> with no atomic relaxation 
were obtained, and the difference of those energies with 
respect to the energy of a system where the tubes and 
surface are 10 A apart are reported (this is called the 
energy gain and it is not an adsorption energy as the 
structures are not fully relaxed). The results are shown 
in Fig. |S1 where the lower curves (red squares) refer to 
SWNTs in the perpendicular configuration. The upper 
curves (blue triangles) refer to the SWNTs in the parallel 
configuration. For our chiral tubes, the range of angles 
is given by (0, 27r/GCD(m, n)), where GCD(m, n) is the 
greatest common divisor of m and n. This range gives 
the maximum angular freedom that the nanotubes will 
have before the position of their carbon atoms becomes 
periodic. Notice the almost complete independence on 
(f) for this energy for all semiconducting tubes, due to 
their weaker bonding as compared with metallic nan- 
otubes. The (6,2) nanotube with the smallest studied 
diameter shows the most pronounced oscillations. This 
is due to its close proximity to the Si(lOO) surface that 
can greatly modify its energy as carbon and silicon atoms 
are brought closer to each other as a result of the tube's 
rotation about its own axis. The distinctive trend for 
the metallic tube indicates an optimal configuration in 
which an extra bond (for a total of eight bonds) might 
form when full relaxation is performed. All the previous 
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FIG. 7: (Color online) Wavefunctions depicting F-point states right below (highest occupied electronic state) and above (lowest 
unoccupied electronic state) the Fermi energy. The hydrogen bottom layer is not shown. 



results shown in this paper were obtained at ^ = 0.0. De- 
spite the fact that the metalhc tube might not be in opti- 
mal angular configuration, our results are the first to date 
to show distinctively different trends for semiconducting 
and metallic nanotubes. In fact, once the metallic tube 
is in its optimal configuration, the trends provided here 
will only be accentuated. In the previously studied cases 
(Refs. 0,13,0), the highly symmetric atomic arrangement 
of the achiral (6,6) SWNT with respect to the underlying 
surface silicon atoms results in optimal angular configu- 
rations in which D is a function of (f>. We consider that in 
a more general setting, chirality precludes such a highly 
symmetric configuration from occurring: the relative po- 
sition of carbon atoms closest to the surface with respect 
to surface silicon atoms becomes a complicated function 
of the chiral angle. In Fig. |Hlwe observe a small energy 
dependence per unit length (of the order of 100 meV/A) 
against the tube's angle of rotation cj). This justifies the 
angular sweeping at fixed D for our chiral SWNTs. 

It is also apparent from Fig. |S1 that the metallic nan- 
otube displays an energy dependence on the angle of rota- 
tion twice as large as the semiconducting SWNTs. This is 
consistent with the fact that metallic tubes stay closer to 
the surface in comparison to semiconducting nanotubes 
of the same diameter. This could also be influenced by 
the different surface reconstruction we employed. 

The prominent difference in the adsorption energies 
presented in Table IVlll suggests a mechanism for the sep- 
aration between semiconducting and metallic nanotubes 
in solid phase over this diameter range, provided they 
are always adsorbed in configurations with the greatest 



energy gain. Experimental techniques for nanotube sepa- 
ration to date rely on the chemistry and dipole moments 
of samples in solution.^^^^ 

IV. CONCLUSIONS 

We have studied semiconducting SWNTs adsorbed on 
the Si(lOO) surface. We used the results from the (9,3) 
SWNT, which compare well to previous published re- 
sults for metallic SWNTs ^3, ^ as a test bed for our 
choice of exchange-correlation potential and surface re- 
construction. Semiconducting SWNTs of the diameter 
range studied are placed at an almost constant distance 
to the surface. Those tubes are 0.5 A closer to the sur- 
face when they are above and aligned with the trench 
between adjacent dimer rows, in comparison with any 
other configuration in which the trench and the tube axis 
do not align. We found a weak angular dependence on 
the adsorption energy of the system, and we believe this 
dependence will be further lowered as the length of the 
nanotube's unit cell increases, due to the loss of sym- 
metry in the relative positions of the carbon atoms in 
the SWNT closest to the Si(lOO) surface atoms. For the 
(9,3) SWNT, we found in agreement with Refs.HHthat 
it remains metallic in cither configuration, but the weak 
angular dependence tells us it will remain at an almost 
fixed height D, independent of its angular orientation. 
We also found a smaller adsorption energy for semicon- 
ducting SWNTs in comparison to metallic SWNTs of 
similar diameter, consistent with weaker bonding. 
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FIG. 8: (Color online) Energy gain vs. axial rotation, prior to 
the relaxation cycle and at fixed height. Red squares show re- 
sults in the parallel configuration, while the blue triangles cor- 
respond to the perpendicular configuration. Q = GCD(m,n). 
The results shown here help identify the best angular con- 
figuration for chiral nanotubes, which can not be known a 
priori. The largest adsorption energies occur for tubes paral- 
lel to the trench. The energy dependence on angle of rotation 
is more marked as the nanotube diameter is decreased, since 
is brought in closer proximity with the surface. The metallic 
nanotube shows the most pronounced energy dependence. 



hybrid systems can be extended, located predominately 
over either the siHcon slab or the SWNT, or be localized 
at the interface between the SWNT and the Si(lOO) form- 
ing a one-dimensional conduction channel. In contrast, 
the LUES tend to be more localized towards the surface 
when semiconducting nanotubes are brought into prox- 
imity. Currently, experimental results on this system are 
starting to emerge. J£J We hope that the results provided 
in this paper motivate further experimental work in the 
area, as the properties described in here might be useful 
for electronic as well as opto-electronic applications. 

Systematic studies as the one presented here are truly 
necessary in order to assess robust properties of this hy- 
brid system against properties that might appear for a 
given choice of chiral angles and geometrical configura- 
tion. In our case, the gap reduction appears for both 
configurations, and the bonding lengths are found to 
be consistently larger for semiconducting nanotubes in 
comparison with metallic ones. This implies a reduced 
absorption energy for semiconducting nanotubes on this 
surface. In turn, this might serve to mechanically attach 
metallic nanotubes and release semiconducting ones by 
suitable heating under dry conditions. |2§. [291 
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